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I .  INTRODUCTION 


This  report  describes  and  lists  a  computer  program  designed  to  handle 
pulse  propagation  problems  when  the  propagating  channel  is  the  earth- 
ionosphere  waveguide  and  is  intended  for  use  in  the  elf/vlf  bands.  The  report 
is  an  extension  of  earlier  work  (ref  1),  which  was  restricted  to  pulse 
propagation  of  the  vertical  electric  field  generated  by  a  ground-based 
vertical  eletric  dipole  source.  The  present  extension  makes  allowance  for 
calculating,  at  any  height  within  the  guide,  all  electric  field  components 
generated  by  electric  dipole  sources  of  arbitrary  orientation  and  elevation. 
Inputs  are  mode  data  (i.e.,  eigenangles  and  excitation  factors)  as  a  function 
of  frequency  as  determined,  for  example,  by  the  waveguide  program  of 
reference  2.  The  mode  data  are  interpolated  by  using  cubic  splines.  The 
pulse  shape  integral  (which  is  a  Fourier  transform)  is  calculated  by 
employing,  at  the  user's  option,  either  a  fast  Fourier  transform  technique  or 
a  Filon  technique.  Normally,  the  fast  Fourier  transform  (FFT)  technique  is 
used.  Advantages  and  disadvantages  of  the  FFT  have  been  discussed  by  Seyler, 
Block,  and  Flynn  (ref  3).  Its  major  advantage  is  a  savings  in  computational 
time,  whereas  a  disadvantage  may  be  that,  strictly,  only  periodic  pulse  trains 
may  be  analyzed.  Thus,  when  a  nonperiodic  pulse  is  considered,  it  must  be 
treated  as  a  periodic  pulse  train  with  period  much  greater  than  the  pulse 
width  in  order  to  obtain  adequate  resolution.  Another  disadvantage  of  the  FFT 
is  that  there  is  no  measure  of  the  accuracy  of  the  integral  evaluation.  For 
this  reason,  a  second  integration  routine  based  on  the  Filon  method  (ref  A)  is 
included.  The  method  is  more  direct  but  much  slower  than  the  FFT.  In  addi¬ 
tion  to  the  integral  evaluation,  output  of  the  Filon  integration  contains  an 
indication  of  the  accuracy  of  the  evaluation,  and  this  is  perhaps  most  useful 
for  purposes  of  checking  the  FFT. 

At  present,  the  program  is  designed  to  handle  only  laterally  homogeneous 
waveguides.  It  is  likely  that  the  subroutine  "CHANEL"  could  be  extended  to 
allow  for  lateral  inhomogeneity  of  the  guide  via  WKB  or  mode  conversion 
methods.  Whereas  the  program  of  reference  1  was  developed  primarily  as  a  tool 
for  calculating  slow-tail  atmospheric  wave  forms  (i.e.,  wave  shapes  in  the  elf 
band  generated  by  lightning  discharges),  the  present  program  is  intended  more 
as  an  aid  to  elf/vlf  system  designers.  Thus,  in  addition  to  the  slow-tail 


waveform  capability,  the  program  allows  for  the  study  of  the  distortion  of 
square-wave  and  Gaussian  pulse  envelopes  (as  well  as  sequences  of  such  pulses) 
and  for  the  analysis  of  spread  spectrum  systems  (ref  5,  6).  However,  it  is 
stressed  that  alternative  input  waveforms  can  be  accommodated  by 
straightforward  modification  of  the  subroutine  XMTR. 

The  mathematical  problem  at  hand  simply  reduces  to  the  calculation  of  a 
Fourier  integral  for  which  the  integrand  is  made  up  of  a  transmitter  spectrum, 
receiver  spectrum,  and  channel  spectrum,  each  of  which  is  discussed  in  the 
following  section.  In  section  III,  the  program  input  is  described  and  in 
section  IV  the  program  structure  is  outlined.  Section  V  contains  output 
description  and  sample  results.  The  appendix  contains  a  program  listing. 


II.  SOURCE,  RECEIVER,  AMD  CHANNEL  MODELS 


OUTPUT  WAVEFORM 

In  the  following,  x,  y,  z  is  a  Cartesian  coordinate  system  with  x  -  z 
being  the  plane  of  propagation  and  z  directed  into  the  ionosphere  with  the 
ground  at  z  =  0.  In  terms  of  the  source,  receiver,  and  channel  spectrums,  the 
output  waveform,  G(x,  zR,  t;  z^) ,  at  a  great  circle  range  x,  altitude  zR,  and 
time  t  generated  by  a  source  at  x  =  y  =  0,  z  =  z^  may  be  written  as  (i  =  V"l) 


G(x,  zr ,  t;  zT) 


-tf 

J  -0 


iiot 


—  |  s (u>)r(iu)h(u>,  x,  zR;  zT)e  dw 


s(w)r(u»)h(u>,  x,  zR;  zT)elwtdu» 


i: 


=  2  Re  |  S(f )R(f )H(f ,  x,  zR;  zT)e2"iftdf 


where 


(1) 


S(f)  =  s(u»)  =  s(2nf)  source  (current  moment)  spectrum  (2) 

R(f)  =  r(u>)  =  r(2nf)  receiver  spectrum  (3) 

H(f,  x,  zR;  zT)  =  h(u),  x,  zR;  zT) 

=  h(2nf,  x,  zR;  z^,)  channel  spectrum.  (4) 

The  second  and  third  equalities  in  equation  1  follow  from  the  requirement  that 
G  be  a  real  quantity,  so  that 


S(f )  =  S*(-f),  R(f )  =  R*(-f),  H(f ,  x,  zR;  z ?) 

=  H*(-f,  x,  zR ;  zT)  (5) 

where  the  asterisk  denotes  the  complex  conjugate.  G  can  represent  any  of  the 
electric  field  components,  Ex>  E^,  E z,  generated,  as  mentioned  above,  by  an 
arbitrarily  oriented  electric  dipole  at  x  =  y  =  0,  z  =  z^,.  The  receiver, 
source,  and  channel  functions  are  described  below. 


RECEIVER 


RECVR  is  a  subroutine  that  can  be  easily  modified  or  replaced  to 
accommodate  the  individual  user’s  needs.  The  particular  RECVR  subroutine 
contained  in  the  program  listing  in  the  appendix  assumes  a  receiver  function 
of  the  form 


R(f)  = 


R(f)  = 


(1 


+  i(f  +  f2)/f3 


t  0 


(i  +  if/f3)'Q  ;  f2  =  o 


(6) 


where  the  frequencies  f^,  f^,  f^  and  the  exponents  P  and  Q  are  read  into  the 
program  via  namelist.  This  receiver  function  allows  for  modeling  receivers 
representative  of  spread  spectrum  systems  (ref  5)  as  well  as  those  used  in 
slow  wave  tail  studies  of  atmospherics.  Observe  that  the  function  satisfies 
the  condition  specified  in  equation  5. 

TRANSMITTER 


TRXMTR  is  a  subroutine  that,  too,  can  be  readily  altered  to  meet  specific 
needs  of  the  user.  In  the  present  program,  four  source  functions  are 
available.  They  are  called  by  setting  IFLGTR  =  1,  2,  3,  or  4.  The  source 
functions  are: 

i)  IFLGTR  =  1 

s(u>)  =  u(ui)v(iu)  (7) 

where 


u(u>) 


iu)T/2 

e _  1 

2  (u>  -  w) 


(■ 


i(u>  -u>)T 
o 


(u>  +  U)) 


(■ 


-i(u>  +w)T 
o 


(8) 


I, 


■iu>n(T+6t) 


im(N  +  1)(T  +  6t)/2)  sin(u»(N  +  1)(T  +  6t)/2 


exp(-iu>(T  +  6t)/2J  sin|u>(T  +  6t)/2) 

The  spectrum  s(u>)  corresponds  to  the  time  function 


where 


^(t)  =  sin|u>o(t  -  n(T  +  6t)  +  T/2)}  :  n(T  +  fit)  -  T/2  <  t  <  n(T  +  fit)  +  T/2 
=  0  :  otherwise. 

Equation  10  represents  a  sine  wave  modulated  by  a  series  of  N  +  1  square 
wave  envelopes  as  indicated  below. 


T/2  T/2  +  fit 


3/2  T  +  fit 


N(T+«t)  -T/2  N(T+Jt)  +T/2 


If  equation  1  is  used  to  define 


“'I 


a  +  ib  =  S(f )R(f )H(f ,  x,  zR;  zT)e2niftdt  , 


the  envelope  +  b2  normalized  to  unity  is  the  plotted  output  for  IFLGTR 

1.  Printed  output  in  units  of  dB  above  a  pV/m  per  kW  is  also  available.  T 


»JA 


normalization  assumes  that  equation  10  is  multiplied  by  a  current  moment 
corresponding  to  a  CW  power  output  of  1  kW  at  fQ  when  placed  vertically  over  a 
perfectly  conducting  half  plane. 

ii)IFLGTR  =  2 


s(uj)  =  Wj^VjCu))  +  w  (u»)v  (u») 


where 


WjCui)  =  ^  exp^-(u»Q  -  w)2T2/4^ 
w2(ui)  =  exp^-(u)Q  +  w)2T2/4^ 


Zi(u>  -u»)n6t 
e  ° 

n=o 

expfi(u)  -  ui)(N  +  l)6t/2)  sinf(u>  -  ui)(N  +  l)6t/2) 


exp^i(uiQ  -  u>)6t/2^ 

sin^(wo  -  u>) fit / 2) 

t  -i(u)  +w)n6t 

e 

f 

(-i(u)Q  +  u>)(N  +  l)6t/2) 

sin((iuo  +  u»)(N  +  l)6t/2) 

The  spectrum  in  this  instance  corresponds  to  the  time  functions 
N 

gj(t)  q,„(t) 


*ln(t> 


where 


q.  (t)  =  exp(-(t  -  n6t)2/T2)sin(wot) 


In  terms  of  the  definitions  for  a  and  b  given  in  equation  11,  the  envelope 

2~Jaz  +  b2,  normalized  to  unity,  is  the  plotted  output  for  IFLGTR  =  2.  Printed 
output  in  units  of  dB  above  a  pV/m  per  kW  is  also  available.  This 
normalization  assumes  that  equation  17  is  multiplied  by  a  current  moment 
corresponding  to  a  cw  power  output  of  1  kW  at  fQ  when  placed  vertically  over  a 
perfectly  conducting  half  plane. 


iii) IFLGTR  =  3 


Rothmuller  (ref  5)  and  Kelly  et  al.  (ref  6)  have  investigated  the  effect 
tnat  the  earth-ionosphere  waveguide  has  on  one  type  of  vlf  communication 
system.  The  system  studied  was  characterized  by  a  differential  phase-encoded 
signal  waveform  composed  by  frequency  shift  keying  (FSK)  a  carrier  with  a 
binary  pseudorandom  or  pseudonoise  (PN)  sequence  of  pulses  or  chips.  The  FSK 
modulation  index  is  0.5,  which  is  designated  as  minimum  shift  keying  (MSK) . 
For  more  detail  concerning  the  basic  waveform  and  terminology,  the  interested 
reader  should  see  references  5  and  6.  Here  we  note  only  that  the  PN  sequence 
has  a  power  spectrum  (or  source  spectrum  for  our  purposes)  given  by 


8  cos2((f  -  f0)2n/fc) 


fc  M1  -  16(f  -  V  /f 


cos2  ((f  +  fQ)2Ti/fc) 


where 


f  =  carrier  frequency 
f  =  chip  frequency 

and  K  is  a  constant  determined  rather  arbitrarily  from  the  relation 


I 


00 

S(f )df  =  K 

-00 


(20) 


K  is  chosen  so  that,  when  used  in  conjunction  with  the  channel  function  given 
subsequently,  it  would  correspond  to  a  vertical  electric  dipole  current  moment 
at  frequency  f^,  which  would  radiate  1  kW  of  power  when  placed  over  a 
perfectly  conducting  plane.  This  normalization  gives 


2.386  x  ip8 
f 

o 


A/m  . 


(21) 


with  f  in  Hz. 
o 


Output  of  the  correlation  receiver  corresponding  to  the  delay  time  T  is 


c 


2Re^l  S(f )R(f )H(f ,  x,  zR;zT)e27llfZ  }=  2  Re(a  +  ib) 
o 


[ifz  |  _ 


=  2  V(a')2  +  (b')2  cos(2nf  T  +  $)  (22) 

o 


where 


i'  =  a  cos  2nf  t  +  b  sin  2nf  x 
o  o 


b'  =  -a  sin  2/if  x  +  b  cos  2nf  X 


(23) 


tan<|>  =  b'/a' 


The  envelope  2V(a')2  +  (b ’ ) 2  expressed  in  dB  above  1  uV/m/kW  radiated 
(interpreted  in  the  sense  of  equations  20  and  21)  and  the  phase,  <Ji,  of  the 
correlation  vector  as  a  function  of  the  delay  X  is  the  output  corresponding  to 
IFLGTR  =  3. 
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iv)IFLGTR  =  4 

The  principal  motivation  for  the  earlier  work  (ref  1)  was  to  study  the 
shape  of  slow  wave  tails  associated  with  atmospheric  discharges.  This 
capability  is  retained  in  the  present  report.  For  this  purpose,  the 
particular  source  function  contained  in  the  subroutine  TRXMTR  is  Williams' 
(ref  7)  mean  source  description  for  a  lightning  discharge,  which  is  given  by 


Idl(ui)  =  v 


y'  _A 

L*  ll,n  +  J 


exp  -  T  (v  +  ju>) 

_ L _ P-n. _ J 

1  +  Xv(Yn  + 


where 


=  16.8  x  io3  A 
=  15.35  x  io3  A 

3 

=  10  A 

=  0.45  x  io3  A 
=  43  (Js 

=  3.5  x  io7  m/s 


Yj  =  5.88  x  io5  s”1 

Y2  =  3.03  x  io4  s"1 

Y3  =  2.0  x  io3  s'1 

Y^  =  1.47  x  io2  s"1 

T  =  180  ps 
v 


The  A!s,  y’s»  X  ,  T  ,  and  v  are  contained  in  DATA  statements  in  TRXMTR.  The 
i  t  p  v  o 

units  of  amperes  for  the  expressions  and  m/s  for  v  ,  coupled  with  the 
channel  defined  in  the  following  subsection,  yield  a  plotted  wave  form  in 
units  of  pV/m.  The  printed  output  is  in  units  of  dB/pV/m. 


CHANNEL-EXCITATION  FACTORS  AND  HEIGHT  GAINS 


Summarized  in  this  subsection  are  modal  excitation  and  height  gain 
formulas  required  as  input  for  the  mode  sum  evaluations,  which  allow  for 
arbitrary  elevation  (within  the  guide)  and  orientation  of  the  transmitter  and 
receiver.  The  formulas  have  been  given  earlier  (ref  8)  and  are  included  here 
for  completeness.  The  excitation  factor  formulas  are  given  in  the  3*3 
matrix  below.  The  column  headings  apply  to  excitation  of  the  electric  field 


9 


components  E  ,  E  and  E  and  the  row  headings  apply  to  excitation  by  a 
z  y  x 

vertical  electric  dipole  (A  ),  horizontal  broadside  electric  dipole  (AB),  and 

V  D 

a  horizontal  end-on  electric  dipole  (Ag).  The  direction  of  z  is  taken 
positive  into  the  ionosphere  with  x  -  z  being  the  plane  of  propagation  and  y 
normal  to  the  plane  of  propagation. 


FIELD 

COMPONENT 


|XCITER 


s2t. 


-ST„ 


ST, 


^E 


-c.t  T 
"34 

-ST, 


-T  T 
3  4 


(26) 


where 


T,  = 


(1  +  „Rtl)2(l  -  ARA  ^RJS 

§1  r  n 

30  "  "  11 


1/2 


(27) 


T2  = 


(1  +  1R1)2(1  -  „R„  „R„)S1/2 

—  R  D 
00  u22 


(28) 


T3  = 


(1  +  „R„)  ( 1  + 

30  12 


1/2 


(29) 


T* =  ^ 


(30) 


The  R  and  R's  represent,  repectively,  elements  of  the  reflection  matrix 
looking  into  the  ionosphere  and  toward  the  ground  from  the  same  level  d  within 
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the  guide.  The  first  subscript  refers  to  the  polarization  of  the  incident 
wave,  while  the  second  applies  to  the  polarization  of  the  reflected  wave.  S 
is  the  sine  of  an  eigenangle,  0,  and  3F/30  is  the  derivative  of  the  modal 
function  which  is  evaluated  at  an  eigenangle.  The  T  and  0  expressions  are 
input  from  the  waveguide  program  of  reference  2. 


The  excitation  factors  must  be  supplemented  with  definitions  of  height 
gains.  These,  along  with  the  definitions  of  the  expressions,  are 


f„(z)  =  exp^Z-  a--^(F1h1(q)  +  F2h2(q)) 


fx(z)  =  F3h1(q)  +  F4h2(q) 


8(2)  =  ik  h  f-(z) 


D„  =  f„2(d)  D12  =  f„(d)fx(d)  D22  =  fx2(d) 

r,  -  -  |m«0>  -  *  r*  (rf'Xv  -  ^'W} 

n  2  /  vl/3/  \l/2 

f2  - «.(«.)  -  i  r*(r)  (Vs7  M,.) 

F,  =  -  |  hJCV  -  ‘(r)1/3(»g2  -  S2)1/2b2Cq0)| 


F«  -  hi(q0)  -  i(r)1/3(v  '  s2)1/2h2(,o) 


*  fen 


C2  "  I  (h  •  2>) 


i 


H.(q>  =  h'(q)  ♦  h.(q) 


j  =  1,2 


(40) 


n2  =  1  -  -(h  -  z) 
a 


(41) 


where 


cosine  of  the  angle  of  incidence  at  height  h 

free  space  wave  number 

dielectric  constant  of  the  ground 

ground  conductivity 

circular  radio  frequency 

earth's  radius. 


(42) 


I 


The  functions  h^  and  h2  are  modified  Hankel  functions  of  order  1/3  (which 

are  linearily  related  to  Airy  functions),  as  defined  by  the  computation 

Laboratory  at  Cambridge,  Massachusetts  (ref  9).  The  primes  on  these 

quantities  denote  derivatives  with  respect  to  the  argument.  Equation  41  is 

the  modified  refractive  index  which  equals  unity  at  height,  h.  The  subscript, 

2 

o,  which  appears  on  n  in  equations  35  and  36,  signifies  that  equation  41  is  to 

be  evaluated  for  z  =  0.  Similarly,  the  subscript  o  that  appears  on  q  in 

equations  35  through  38  signify  that  equation  39  is  to  be  evaluated  for  z  =  0. 

It  should  be  pointed  out  that  f„  is  the  height  gain  for  the  vertical  electric 

field  E  ,  fx  the  height  gain  for  the  horizontal  electric  field  component  (E  ) 
z  y 

normal  to  the  plane  of  propagation,  and  g  the  height  gain  for  the  horizontal 
electric  field  component  (E^) ,  which  is  in  the  plane  of  propagation. 


Because  the  imaginary  part  of  the  eigenangle  in  absolute  value  can  become 
quite  large  when  operating  in  the  ELF  range,  it  proves  necessary  to  avoid 
overflow,  and  use  of  the  flat  earth  analogues  of  equations  31  through  33  is 
justified.  That  is,  to  replace  the  height  gains  by 


f„(z)  =  exp(ikCz)  +  „RtI  exp(-ikCz  +  2ikCd) 


(43) 


fA(z)  =  exp(ikCz)  +  ARA  exp(-ikCz  +  2ikCd)  (44) 

g  =  C  jexp(ikCz)  -  „R„  exp(-ikCz  +  2ikCd)J  .  (45) 

When  the  absolute  value  of  the  imaginary  part  of  the  eigenangle  exceeds  10°, 
the  height  gain  functions  will  be  computed  by  equations  43,  44,  and  45. 

The  punched  output  of  the  waveguide  program  of  reference  2  is  transformed 
to  correspond  to  d  =  0,  independent  of  the  actual  d  used  in  the  waveguide  run. 
Therefore,  in  the  present  program  d  in  the  above  formulas  is  set  to  zero. 

CHANNEL-MODE  SUM 


In  terms  of  the  excitation  factors  and  height  gains  defined  in  the 
previous  section,  the  mode  sum  for  the  laterally  homogeneous  guide  may  be 
written  as  follows 


E.(x)  =  - V 

J  jsin(x/2)J  ^ 


|a"  cos(Y)f”(zT)  +  sin(v)sin(<|»)fA(zT) 


+  *E  sin(T)cos(<|i)gn(zT)|  f”(zR) 


-ik(S  -  l)x 
n 


(46) 


The  transmitter  coordinates  are  (0,  0,  z^)  and  the  receiver  coordinates  are 
(x,  0,  2^).  The  mode  index  is  n  and  the  index  j  takes  on  three  values 
corresponding  to  the  electric  field  component  measured  at  the  receiver. 


j  =  1  ■*  z  component  ■*  f  ^  =  f„ 

j  =  2  ■+  y  component  ■*  f  =  fA 

j  =  3  •*  x  component  -*  f^  =  g  • 


M  is  the  dipole  moment  in  A/m  and  Q  is 
Q  =  9.023  x  10_8(if)3/2 


(47) 


with  f  in  Hz.  With  this  value  for  Q,  the  field  is  in  units  of  pV/m.  If  M  is 
chosen  to  correspond  to  1  kW  radiated  power  when  the 


III.  DESCRIPTION  OF  INPUT 


All  input  to  the  pulse  shape  program  is  read  in  via  the  card  reader.  A 
listing  of  a  sample  input  showing  the  data  deck  setup  is  given  on  pages  20 
through  23. 

There  are  two  parts  to  the  input.  The  first  part  is  read  in  by  means  of 
an  ASCII  FORTRAN  name  list  input  format.  The  name  list  read  in  is  initiated 
with  a  control  card  NAME  in  columns  1  through  4.  The  following  variables  and 
arrays  may  be  specified  in  the  name  list  input. 

NFFT  -  2**NFFT  is  the  number  of  integration  intervals  in  the  range 

(FREQU-FREQL)  when  using  the  FFT  for  the  Fourier  evaluations. 
Default  value  is  11,  which  is  the  maximum  allowed  without 
changing  dimensioning. 

FREQU  -  upper  frequency  of  integration  in  kilohertz.  Default  value  is 
100  kHz. 

FREQL  -  lower  frequency  of  integration  in  kilohertz.  Default  value  is 
0.0  kHz. 

INTPRT  -  flag  to  control  print  interval  for  transmitter,  receiver, 

channel,  and  product  spectra  as  a  function  of  frequency.  For 
example,  the  first  NPRNT  (see  below)  spectra  are  printed  and 
thereafter  only  those  spectra  for  which  the  frequency  index 
modulo  INTPRT  equals  zero  will  be  printed.  Default  value  is 
100. 

NPRNT  -  flag  to  control  print  interval  for  transmitter,  receiver, 

channel  and  product  spectra  as  a  function  of  frequency.  The 
first  NPRNT  spectra  are  printed.  Default  value  is  40. 

TAUMAX  -  controls  the  latest  time  in  seconds  plotted  on  the  output 

waveform  curve.  Also  controls  the  maximum  time  for  which  the 


Fourier  integral  is  calculated  when  using  the  Filon  routine. 
Default  value  is  0.002  s. 


TAUO 

NUMTAU 

FREQO 

PULSEW 

FREQ1,  FREQ2 , 
FREQ3 

P,  Q 

RHOMIN 

RHOMAX 

DELRHO 


controls  the  earliest  time  in  seconds  plotted  on  the  output 
waveform  curve.  Also  controls  the  minimum  time  for  which  the 
Fourier  integral  is  calculated  when  using  the  Filon  routine. 
Default  value  is  -0.001  s. 

number  of  taus  between  and  inclusive  of  TAUO  and  TAUMAX  for 
which  Filon  evaluations  are  to  be  made.  Default  value  is  41. 

carrier  frequency  in  kHz.  Default  value  is  23  kHz. 

square-wave  pulse  width  or  Gaussian  e  1  half  width  in  ps. 
Default  value  is  600  ps. 


frequencies  in  kHz  appearing  in  the  receiver  function  given  by 
equation  6.  Default  values  are  FREQ1  =  0.01  kHz,  FREQ2  = 

0.0  kHz,  FREQ3  =  2.5  kHz. 

exponents  appearing  in  the  receiver  function  given  by  equation 
6. Default  values  are  P  =  0.0,  Q  =  2.0. 

minimum  range  in  km  for  which  mode  sum  and  pulse  shape  or 
correlation  vector  is  to  be  evaluated.  Default  value  is 
1000  km. 

maximum  range  in  km  for  which  mode  sum  and  pulse  shape  or 
correlation  vector  is  to  be  evaluated.  Default  value  is 
1000  km. 


incremental  ranges  in  km  between  RHOMIN  and  RHOMAX  for  which 
mode  sum  and  pulse  shape  or  correlation  vector  is  to  be 
evaluated.  Default  value  is  1000  km. 
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TALT 


RALT 


INCL 


THETA 


I  COMP 


IFLGTR 


INTFLG 


CHIPFR 


NUMPLS 


transmitter  altitude  in  km.  Default  value  is  0,  which 
corresponds  to  ground-based  transmitter. 

-  receiver  altitude  in  km.  Default  value  is  0,  which  corresponds 
to  ground-based  receiver. 

-  inclination  of  transmitter  from  positive  z  axis  in  degrees 
(angle  y  in  equation  44).  Default  value  is  0°,  which 
corresponds  to  a  vertical  antenna. 

-  azimuth  of  transmitter  measured  counterclockwise  from  x  axis  in 
degrees.  Default  value  is  0° ,  corresponding  to  an  end  on 
launch. 

-  singles  out  electric  field  component  calculated.  ICOMP  =  1 

gives  th^  vertical  or  E  field.  ICOMP  =  2  gives  the  E  field 

z  y 

and  ICOMP  =  3  gives  the  Ex  component.  Default  value  is  1. 

-  selects  input  wave  form.  IFLGTR  =  1  gives  sequence  of  square 
wave  pulses.  IFLGTR  =  2  gives  sequence  of  Gaussian  pulses. 
IFLGTR  =  3  gives  waveform  composed  by  frequency  shift  keying 
(FSK)  a  carrier  with  a  binary  pseudorandom  or  pseudonoise  (PN) 
sequence  of  pulses  or  chips.  The  FSK  modulation  index  is  0.5. 
IFLGTR  =  4  gives  Williams'  source  for  the  slow  wave  tail 
calculation.  Default  value  is  1. 

-  selects  integration  scheme.  INTFLG  =  0  invokes  FFT  algorithm 
and  is  the  normal  operating  mode.  INTFLG  f  0  invokes  the  Filon 
method.  Default  value  is  0. 

-  chip  frequency  (f^  in  eq  19)  in  kHz  to  be  used  with  IFLGTR  =  3. 
Default  value  is  1.  kHz. 

-  number  of  pulses  in  the  square  wave  or  Gaussian  sequence. 
Default  value  is  1. 
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PULSED  -  time  delay  in  us  between  square  wave  or  Gaussian  pulses. 

Default  value  is  600.  ps. 

IPLOT  -  flag  controlling  plots  of  transmitter  and  receiver  spectra. 

IPLOT  =  0  gives  no  plots.  IPLOT  ^  0  gives  both  spectrum  plots. 
Default  value  is  0. 

IPL0T1  -  flag  controlling  plots  of  channel  and  product  spectra.  IPL0T1 
=  0  gives  no  plots.  IPL0T1  f  0  gives  both  spectrum  plots. 
Default  value  is  0. 

Following  the  namelist  input,  the  second  part  of  the  input,  consisting  of 
run  identification  and  mode  constant  cards,  is  read  in.  The  read  in  is 
initiated  by  a  control  card  with  DATA  in  columns  1  through  U.  Run 
identification  is  read  in  by  using  a  20A4  format.  Mode  data  appearing  on 
pages  20  through  23  consist  of: 

R  -  punched  output  from  waveguide  program,  which  is  not  used  in  present 
program. 

F  -  frequency  in  kHz. 

A  -  azimuth  in  degrees  relative  to  geomagnetic  north,  for  which 
waveguide  program  was  run. 

C  -  codip  in  degrees  of  geomagnetic  field,  for  which  waveguide  program 
was  run. 

2 

M  -  strength  of  geomagnetic  field  in  Weber/m  ,  for  which  waveguide 
program  was  run. 

S  -  ground  conductivity  in  Siemens/m. 

E  -  dielectric  constant  of  the  ground. 


T  -  punched  output  from  waveguide  program,  which  is  not  used  in  present 
program. 

The  card  containing  the  above  information  is  followed  by  two  cards  for  each 
mode  at  frequency  F.  The  first  card  contains  the  index  1  and  the  following 
mode  input  data : 

TR1  -  real  part  of  the  eigenangle  in  degrees. 

Til  -  imaginary  part  of  the  eigenangle  in  degrees. 

TMP1  -  excitation  factor  coefficient  T^  given  in  equation  27. 

TMP2  -  excitation  factor  coefficient  T2  given  in  equation  28. 

The  second  card  contains  the  index  2  and  the  following  mode  input  data: 

TR1  -  real  part  of  the  eigenangle  in  degrees  (repeat  of  input  of 

information  on  the  first  card). 

Til  -  imaginary  part  of  the  eigenangle  in  degrees  (repeat  input  of 

information  on  the  first  card). 

TMP3  -  excitation  factor  coefficient  T^  given  in  equation  29. 

TMP4  -  excitation  factor  coefficient  T^  given  in  equation  30. 

After  all  data  are  read  in  for  a  given  frequency,  mode  data  for  the  next 
frequency  are  read  in,  etc.  After  data  for  all  of  the  frequencies  are  read 
in,  transfer  of  the  program  to  its  execution  phase  is  initiated  by  a  control 
card  with  STAR  in  columns  1  through  4.  Upon  completion  of  the  program  for  the 
first  complete  set  of  data,  a  new  complete  set  of  data  can  be  read  in, 
processed,  and  executed,  etc. 


SAMPLE  INPUT 


3 

1 

3 

9 

ki 

3 

3 

I NT PR T  «  : 

4 

TAUMAX  »  : 

5 

IFLGTR  « 

6 

NPRNT  «  21 

7 

TAUO  =  -1 

8 

NUMPlS  ■  : 

9 

FREQ2 -23 .  i 

10 

I  PLOT  *1  ,- 

11 

I  PLOT  1*1, 

12 

6END 

13 

DATA 

14 

BETA=0 . 5 ,  1 

15 

R 

.000  F 

16 

1 

89.54299 

17 

2 

89.54299 

18 

1 

83.15179 

19 

2 

83.15179 

20 

1 

76.02858 

21 

2 

76.02858 

22 

1 

71  .53966 

23 

2 

71 .53966 

24 

1 

64.67405 

25 

2 

64.67405 

26 

1 

60.43294 

27 

2 

60.43294 

28 

29 

R 

.000  F 

30 

1 

89.87716 

31 

2 

89.87716 

32 

1 

89.42791 

33 

2 

89.42791 

34 

1 

82.12311 

35 

2 

82.12311 

36 

1 

78.92260 

37 

2 

78.92260 

38 

1 

74.36239 

39 

2 

74.36239 

40 

1 

71  .51556 

41 

2 

71  .51556 

42 

1 

66.95303 

43 

2 

66.95303 

44 

1 

64.15047 

45 

2 

64.15047 

46 

1 

59.43416 

47 

2 

59.43416 

48 

1 

56.38749 

49 

2 

56.36749 

SO 

51 

R 

.000  F 

52 

1 

89.95078 

S3 

2 

89.95078 

54 

1 

89.72677 

.  667, P*0. 0,0  =  2 


.331871 


.972492 

.420081 


15.0000  A 
-4.447931 
-4.447931 
-2.627812 
-2.627812- 
-.339621- 
-.339621 
-.345502 
-.345582- 
-.423671- 
-.423671 
-.501782 
-.501782- 
-.501101- 
-.501101 
-.766432 
-.766432- 
-.548001- 
-.548001 
-1 .181612 
-1 .181612- 


96.648  C 
5.93769233- 
1 .71712808- 
1 .40675269- 
4.75598456- 
-1.811 76863 
6.34841966- 
3.57218043- 
-1 .48832975- 
-5.7463161 1- 
5.05197228- 
8.06905213- 
-3.41322277- 

96-848  C 
7.97650921* 
6.98923905- 
1 .81 225575- 
-1.05341859- 
-1.10011261- 
1  .4102391(1- 
2.1 3041 180- 
-2.25597727- 
-2.86681467- 
4.13103759- 
4.20202594- 
-5.38389344- 
-5.60380652- 
9.72961 132- 
7.26349233- 
-1 .14101552- 
-8.34006036- 
1 .91310056- 
1 .03418070- 
-2.09625628- 


26.610 
-004-4 . 
-008  1 . 
-003-4 . 
-010-4. 
-003-5. 
-008  1 . 
-003-9. 
-009-1 . 
-003-5. 
■007  2. 
-003-3. 
-007-3. 

26.610 
-006-1 . 
-008  1 . 
-003-1 . 
-007-3. 
-003-4. 
-007  5. 
-003-1 . 
-007-8. 
-003-3. 
-007  1 . 
-003-3. 
-007-1 . 
-003-3. 
-007  2. 
-003-6. 
-006-3. 
-003-2. 
-006  3. 
-002-1 . 
-006-4. 


M  4.680- 
33237092- 
76637679- 
30542375- 
42039628- 
79340308- 
00232208- 
53919356- 
55049825- 
44617977- 
73988729- 
46847059- 
64204141- 

M  4.680- 
86808831- 
73850939- 
87533820 
1 061 8208- 
10230709- 
33180142- 
65751675- 
10605862- 
58314899- 
35826212- 
25397411- 
75630743- 
21027605- 
50720018- 
124093G6- 
06261 120- 
83699052- 
72124325- 
02234199- 
50653664- 


005  S 
002-2. 
007  9. 
004-1 . 
007  9. 
002-4. 
006  9. 
004-6. 
006  9. 
002-5. 
006  9. 
003-1  . 
006  9. 

005  S 
002-1 . 
007  9. 
003-3. 
007  9. 
002-3. 
007  9. 
003-2. 
007  9. 
002-3. 
006  9. 
003-6. 
006  9. 
002-1  . 
006  9. 
003-1  . 
006  9. 
002-3. 
006  9. 
002-1 , 
006  9. 


4 . 640+000 
66095085- 
76119950- 
30534744- 
76321563- 
57239026- 
76755522- 
43127843- 
77606788- 
96155712- 
79136109- 
51 790046- 
81437340- 


E  81 
013-6 
001-1 
010-1 
001-1 
012-1 
001-1 
010-5 
001-1 
01  1-1 
001-1 
009-6 
001-1 


4 . 640+000  E  81  . 
36770776-012-1  . 
54479575-001-9. 
94927788-011-2. 
54131447-001-9. 
93972420-012-5. 
52227563-001-9. 
50544481-010-4 . 
50611316-001-9. 
07619798-01 1-4. 
47571643-001-9. 
03080086-010  2. 
45343502-001-1 , 
31560195-010-1  . 
4132351 1 -001-1  . 
05415335-009  1. 
39071193-001-1 . 
93695111-010-3. 
34739329-001-1 . 
56078298-009  3. 
33384195-001-1 . 


26.610  M  4.680 
-004-3.08029371 
-008  1 .00341518 
-003-5.19614783 


■005  S  4.640  +  000  E  81  .0  T 
-003-5 .71833862-012-7 . 4570 
■007  9.70444761-001-3.0651 
■003-7.55528420-012-3.4144 


20 


55  2  89.72677  -4.972712-1.15447172-007-1.70629631-007  9.70000498-001-3.0777 

56  1  85.57822  37453 1 -9 .40 1 72227-004-3 . 2 1 498699-00  2-2 . S3  1 603 19-01 2-3 . 6242 

57  2  85.57822  -.374531  1.30888923-007  3.71334362-007  9.66683730-001-3.0583 

58  1  83.08694  -. 240632  2 . 028 1 2662- 003-2 . 06 1 06337- 003- 1 . 1 547 1 073-0 1 0  5 . 1 048 

59  2  83.08694  -.240632-2.1989547-1-007-5.45562628-007  9.64558974-001-3.0479 

60  1  79.12334  -.397391-2.51086711-003-2.78670846-002-2.04451764-011-2.5083 

61  2  79.12334  -.397391  3.51979452-007  9.05249671-007  9.59300555-001-3.0774 

62  1  77.11312  -.275102  3 . 306 1 397 7-003-3 . 1 8349389-003-3 . 0 1 428 1 01 -01 0  2.1585 

63  2  77.11312  -.275102-4.72247425-007-1.10789632-006  9.55758922-001-3.0714 

64  1  73.44476  -.459411-4.59294498-003-2.38388672-002-8.37238039-011-9.7335 

65  2  73.44476  -.459411  7.61633153-007  1.64506844-006  9.47922789-001-3.1771 

66  1  71.63047  -.3768125.62715752-003-5.50687622-003-5.23718499-0109.5978 

67  2  71.63047  -.376812-9.37261717-007-1.88752817-006  9.43363637-001-3.2154 

68  1  67.93610  -.511771-7.17936538-003-1.97459643-002-2.53456083-010-2.5509 

69  267.93610  -.511771  1.47960783-0062.42963196-0069.33020763-001-3.4446 

70  1  66.14730  -.531512  8.40076327-003-9.21373558-003-7.26669051-010  2.5910 

71  2  66.14730  -.531512-1.73095953-006-2.71428090-006  9 . 27492820-001 -3 . 5861 

72  1  62.41710  -.555411-9.07247630-003-1.51738158-002-6.00887652-010-4.8687 

73  2  62.41710  -.555411  2 . 4670602 3-006  2.97335941-006  9.15089674-001-3.9828 

74  1  60.47119  -.741961  1 . 05556 1 1  I  - 002-1  . 39936338-00  2-B . 50278722-0 1 0  5 . 0284 

75  2  60.47119  -.741961-2.82933254-006-3.32145967-006  9.08521205-001-4.3171 

76  1  56.76943  -. 6 1 0202-9 . 80309330-003- 1 . 11 6 1 6928-002- 1 . 1 33836 16-009-7 . 4866 

77  2  56.76943  -.610202  3.50012604-006  3 . 2 1 396806-006  8.94910932-001-4.9218 

78 

79  R  .000  F  25.0000  A  96-848  C  26.610  M  4.680-005  S  4.640+000  E  81.0  T 

80  1  89.95961  -6.296002  4 . 1 4965702- 005-2 . 84423062-004-2 . 65959360-01 2-5 . 6880 

81  2  89.95961  -6.296002  1.78132298-008  2.11544882-008  9.99736615-001-4.8076 

82  1  89.75731  -5.7023G1  6 . 6G626482- 004-2 . 24545391 -003- 1 . 42592036-01 2-4 . 5403 

83  2  89.75731  -5.702361-2.7598071(1-008-5.24136694-008  9.98948187-001-4.6248 

84  1  88.74632  -.740711-1.14792565-003-2.27329135-002-3.50932488-012-3.4240 

85  2  88.74632  -.740711  1.37327154-007  3.05298592-007  9.95857358-001-2.7427 

86  1  86.21623  -.323792  2.80416722-003-3-36569911-003-5.43727906-011  1.7317 

87  2  86.21623  -.323792-2.14213904-007-4.40146156-007  9.94571388-001-2.1050 

88  1  82.17971  -.352921-2.85400430-003-2.34048092-002-1.54511952-011-1.7795 

89  2  82.17971  -.352921  3.04256552-007  6.64677630-007  9.89881411-001-7.2338 

90  1  80.61306  -.23,5442  3.28247802-003-3.42158275-003-1.68198547-010  1.6271 

91  2  80.61306  -.235442-3.93487753-007-8.10396045-007  9.87260617-001  1.4355 

92  1  77.39001  -. 3758G 1 -4 . 55 1 83 1 t 1 -003- 1 . 89235786-002-6 . 26825926-0 1 1 -6 . 8954 

93  2  77.39001  -.375061  6.39205197-007  1.20104490-006  9.79809105-001  4.2323 

94  1  76.07253  -.2931525.18735842-003-5.37725317-003-2.98127936-0106.9079 

95  2  76.07293  -.293152-7.58520216-007-1.32923716-006  9.76186506-001  5.9032 

96  1  72.92168  -.393581-6.79070043-003-1.46260866-002-1.92001999-010-1.8119 

97  2  72.92168  -.393581  1.23595613-006  1.69775085-006  9.65308584-001  9.1862 

98  1  71.71231  -.392682  7.51004531-003-8.79517768-003-3.95414586-010  1.8395 

99  2  71.71231  -.392682-1.39636245-006-1.81699510-006  9.60419364-001  1.0552 

100  1  68.56366  -.399602-8.11186887-003-9.74748202-003-4.66080674-010-3.3417 

101  2  68.56366  -.399602  2.04390165-006  1.87442600-006  9.45522405-001  1.3631 

102  1  67.33788  -.523181  8.86612176-003-1.34191969-002-3.99816361-0103.3996 

103  2  67.33788  -.523181-2.26887723-006-1.97178016-006  9.38771591-001  1.3018 

104  1  64.22144  -.410722-7.88689498-003-5.51680942-003-8.83920132-010-4.6723 

105  2  64.22144  -.410722  2.77423089-006  1.66314987-006  9.19652522-001  1.5219 

106  1  62.85581  -.666111  8.66894360-003-1.77447477-002-3.15086311-0104.7624 

'07  2  62.85581  -.666111-3.06723743-006-1.72537347-006  9.10496563-001  1.2858 

108  1  59.81876  -.444912-7.08291441-003-2.62661430-003-1.40096686-009-5.7675 

109  2  59.81876  -.444912  3.35079508-006  1.28801229-006  8.88060510-001  1.0387 

110  1  58.20465  -.807741  7.90575973-003-2.10266041-002-1.87648097-0105.8927 

111  2  58.20465  -.807741-3.76355356-005-1,29681153-0068.76700148-001  3.7316 
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0000  A 
801582 
.801582 
.172221 
.  172221- 
.077971- 
.077971 
.802442 
.802442- 
.325201- 
.325201 
.256792 
.256792- 
.314231- 
.314231 
.  287462 
.  287462- 
.303281- 
.303281 
.372092 
.  372092- 
.283522- 
.  283522 
.  48414 1 
.484141- 
.269702- 
.  269702 
.600941 
.600941- 
.272232- 
.272232 
.711151 
.71 1151- 
.294362- 
.294362 
.813321 
.813321- 
.339362- 
.339362 
.90661 1 
,90661 1- 

.0000  A 

,  154792 

.  154792 

.510361 

,518361- 

,  242822 

.  242822 

.  394181 

.394181- 

.313841- 

.313841 

.320462 

.320462- 

.272241- 


96.848  C 
9.21416370 
3.33612060 
2.33388033 
6.15653212 
1 .35654240 
1 .75959459 
3.43013523 
2.37264798 
3.52030885 
2.84209779 
3.77914089 
3.57363756 
4.87920851 
5.79850294 
5.28538768 
6.64902906 
6.47681038 


00696990 
98046980 
19382139 
66976604 
65506199 
1914551 1 
70910638 
52105560 
2.01539578 
6.02721 120 
■2.18522109 
4.25863202 
2.21227225 
4.72565473 
2.43150367 
3.26558537 
2.38145702 
3.65956590 
2.6693371 l 
2.4G994873 
2.60870254 
2.72667289 
■2.98982195 


26.61 
-006-4 
-009  4 
-004-5 
-009-1 
-003-9 
-007  2 
-003-7 
-007-3 
-003-2 
-007  5 
-003-4 
-007-6 
-003-1 
•007  9 
-003-6 
-007-1 
■003-1 
-006  1 
-003-9 
-006-1 
-003-5 
-006  t 
-003-1 
-006-1 
-003-2 
-006  7 
-003-1 
-006-7 
-003-2 
-006  3 
003-1 
-006-2 
-003  9 
-006-1 
-003-2 
-006  3 
-003  2 
-006-7 
-003-2 
-006  1 


0  M  4.680-005 
.  06602093-005 
.35358660-009 
.90389849-004 
.59281 199-008 
.89436009-003 
.  22063330-007 
.  66506500  -003 
.33145127-007 
.  00907697-002 
. 59975376-007 
.37141508-003 
.55233706-007 
.5231 1 159-002 
.  35443637-007 
.  02913019-003 
.01968099-006 
.  05080334-002 
.20814860-006 
.46757186-003 
.  27433299-006 
.72501693-003 
. 12109299-006 
.37139540-002 
. 15557695-006 
.30329942-003 
.  55881238-007 
.693231 10-002 
.43441461-007 
.97679551-004 
.  15480715-007 
.  DD52 1 3B0-002 
.  37438300-007 
.93770670-004 
.71861332  007 
.04236524-002 
.42168367- 007 
.02460142-003 
.50150313-007 
.18211 829-002 
.05839401-006 


S  4 . 640+000 
-7.33408398- 
1 .01745796+ 
■4. 19298588- 
1 .01659237+ 
-7.60637224 
1 .01431617  + 
-2.01 165827- 
1 .01360975+ 
-1 .34843656- 
1 .01064682+ 
-9.63985307- 
1 .00918959+ 
-5.40644682- 
1  .00445074  + 
-1  .77346710- 
1  .  00235936  + 
-1 .66328987- 
9.95344214- 
-2. 19153317 
9.92226020- 
-3.90159866- 
9.81998667- 
-1 .82826587- 
9.76069076- 
-6.97587266 
9.61933985- 
-9.76922293- 
9.53540497- 
-1  .05524151- 
9.32590367- 
1  . 70044362- 
9.20637153- 
•1 .46713304- 
8.949001 13- 
1 .19203678- 
8.81666422 
-1  .95284630- 
8.56219366- 
2.76802057- 
0.45654504- 


E  81  . 
013-2. 
000-1  . 
013-2 . 
000-9. 
012-2 
000  1  . 
011  1  . 
000  1  . 
011-1. 
000  4, 
011  1. 
000  5. 
011-5. 
000  9. 
010  5. 
000  1  . 
010-1  . 
001  1  . 
010  1  . 
001  1  . 
010-2 . 
001  2 
010  2. 
001  2. 
010-2. 
001  3. 
011  2 
001  3 
009-2. 
001  4. 
012  2. 
001  4. 
009-2 . 
001  4. 
010  2  . 
001  3. 
009-3. 
001  2. 
010  2. 
001  6. 


96.848  C 
1 .7195161 1- 
6.26939424- 
6.77289499- 
1 .31226660- 
3.86689535- 
7.92866732- 
1 .97066873- 
1 .22145323- 
3.94388248- 
2.85340914- 
4.36410488- 
3.55639354- 
5.0378381 3- 


26.610 
006-6 . 
010  8. 
005-1  . 
009-4  . 
005-1  . 
008  ff. 
003-8 . 
007-1  . 
003-1  . 
007  4  . 
003-6. 
007-5. 
0C3-1  . 


M  4 . 680-005 
65459652-006 
70561 630-01 0 
3731 1 265-004 
56250016-009 
94379015-003 
55532294-008 
00630264-003 
70501769-007 
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21501258-002 


S  4 . 640+000  E  81 . 
■1 .70822676-013-1 . 
1 .02251831+000-1 . 
-1 .28066713-013-7 , 
1 .02154364  +  000-7. 
-7. 12259201-012-4. 
1 .01944663+000  1 . 
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1 .01870005+000  1 . 
■1 .34537301-01 1-1 . 
1 .01656783+000  3 . 
-5.23506863-01 1  1 . 
1 .01554948+000  3. 
-4.91220250-01 1-4. 


0  T 
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9350 
4027 
7334 
7748 
4141 
2818 
3105 
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3522 
3897 
3568 
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3305 
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3472 
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8041 
6241 
6731 
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6940 
8148 
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7023 
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0  T 
1727 
1474 
8820 
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0835 
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7556 
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0957 
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0193 
9114 
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194 
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196 
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39516 
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60. 

39516 
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58. 

75931 

199 
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58. 

75931 

200 

1 

57. 

22259 

201 
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57. 

22259 

202 
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END 

204 

start 

-.272241  5.43721526- 
-.298492  5.24895213- 
-. 298492-6. t 1081887- 
-.244612-5.64404484- 
-.244612  9.33826691- 
-.365431  5.98990329- 
-.365431-1 .0131 1328- 
-.217262-4.96455305- 
-.217262  1 .26447735- 
-.456311  5.33993304- 
-.45631 1-1 .35897561- 
-. 199072-3.60672452- 
-.199072  1.40992418- 
-.548311  4.04662930- 
-.54831 1-1 .52351 126- 
-.  193242-2.58837634- 
-. 193242  1 .45006429- 
-.634701  2.90467023- 
-.634701-1 .58384565- 
-. 199332-1 .73459633- 
-.199332  1.46168267- 
-.715771  1.97376651- 
-.715771-1 .63377436- 
-.217322-^.96977105- 
-.217322  1.48137963- 
-.791751  1.10823596- 
-.791751-1 .69562365- 
-.248472-2.39132431- 
-.248472  1 .52686277- 
-.861241  1.39434616- 
-.861241-1 .78866654- 
-.296182  6.81471960- 
-.296182  1.60517165- 


007  7.28336296-007  1 
003-7 . 24678044-003-1 
007-7.89565838-007  1 
003-7 . 40134291-003-1 
007  8.43385997-007  1 
003-1 .03043913-002-1 
006-0.05114169-007  1 
003-3.48515713-003-3 
006  6.82723261-007  9 
003-1 .34751227-002-8 
006-6.93080402-007  9 
003-1 . 1 2099801-003-5 
006  3.86070006-007  9 
003-1 .54822026-002-3 
006-3 . 53674 1  18-007  9 
003  1.79750459-004-7 
006  5.46746862-008  9 
003-1 .66300139-002  2 
006  3.66019841-008  9 
003  9.74064271-004-1 
006-3. 18594022-007  9 
003-1 . 74074376-002  9 
006  4.80440357-007  9 
004  1 .55366775-003-1 
006-7.67130373-007  9 
003-1 .B0875934-002  1 
006  1 .04387406-006  9 
004  2.07731486-003-1 
006-1 .33272954-006  8 
004-1 .88362205-002  2 
006  1 .76271519-006  6 
004  2.65892252-003-2 
006-2.07563082-006  8 


.01222657+000  6.5597 
.05158253-010  4.0276 
.01078829+000  7.6739 
.44906575-010-9. 1 144 
.00627670+000  1 . 1729 
.21564057-010  9.2064 
.00417314+000  1 .3086 
. 13408421-010-1 .3144 
.90098433-001  1.9012 
.05600166-01 1  1 . 3299 
.94550772-001  2.0513 
.27402795-010-1 .4485 
.85036089-001  2.0270 
.42179783-011  1 .4585 
.79646951-001  2.9281 
.71907724-010-1 .4315 
.66708526-001  3.7370 
.58675477-011  1.4094 
.57122207-001  3.5085 
.04989406-009-1 .3279 
.39793691-001  4.0768 
.421 13323-01 1  1 .2209 
.28504840-001  3.4118 
.36881209-009-1 . 1086 
.10345115-001  3.2940 
.79295710-010  8.2145 
.01  143454-001  2 . 1334 
. 74026334-009-6 . 7406 
.86461414-001  1.5787 
.94850606-010  2.9668 
.80726665-001  2.4958 
.18246249-009  1.6829 
.70145038-001-3.8564 
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IV.  PROGRAM  LAYOUT 


P* 


This  section  describes  only  the  basic  features  of  the  pulse  shape  program 
listed  in  the  appendix.  In  particular,  plot  and  label  routines  PLSPEC, 

BORDER,  SYMBOL,  CURVE,  PLOT,  WOPLOT,  PL0T12,  PLOT3  and  PLLABL,  are  not 
described. 

MAIN  controls  the  program  flow.  Subroutines  it.  the  order  of  their  call 
are  described  below: 

SUBROUTINE  INPUT 

INPUT,  called  from  MAIN,  reads  in  NAMELIST  and  mode  data. 

SUBROUTINE  HTGAIN(IOPT,  FREQ,  SIGMA,  EPSR,  ALPHA,  NRMODE ,  TP,  Z,  HG) 

HTGAIN,  called  from  MAIN,  evaluates  the  modal  height  gain  functions  that 
appear  in  equation  46.  The  arguments  of  HTGAIN  are: 

IOPT  -  option  flag  set  for  1  in  present  program. 

FREQ  -  frequency  (kHz). 

SIGMA  -  ground  conductivity  (s/m). 

EPSR  -  real  dielectric  constant  of  the  ground. 

-4  -1 

ALPHA  -  earth  curvature  constant  (3.14  x  10  km  ). 

NRMODE  -  number  of  modes  at  any  one  of  the  input  frequencies. 

TP  -  complex  ground  eigenangle  in  degrees, 

z  -  altitude  at  which  height  gain  is  evaluated  (km) . 

HG  -  height  gain. 
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SUBROUTINE  MDHNKL(Z,  HI,  H2,  H1PRME,  H2PRME,  THETA,  IDBG) 


This  subroutine  is  called  by  HTGAIN.  It  evalutes  the  modified  Hankel 
functions  of  order  1/3  and  their  derivatives  according  to  the  methods  of 
reference  9.  The  arguments  of  MDHNKL  are: 

Z  -  argument  of  modified  Hankel  functions  of  order  1/3  and  their 

derivatives . 

HI,  H2  -  modified  Hankel  functions  of  order  1/3. 

H1PRME, 

H2PRME  -  derivatives  of  modified  Hankel  functions  of  order  1/3. 

THETA  -  not  used  in  present  program. 

IDBG  -  not  used  in  present  program. 

SUBROUTINE  FUNSPL(MD.  LF) 

FUNSPL  is  called  from  MAIN.  Inputs  to  FUNSPL  are  a  mode  index,  MD,  which 
takes  on  values  1  through  the  maximum  number  of  modes  read  in  for  a  given 
frequency;  and  the  index,  LF,  for  the  quantity  that  is  to  be  approximated  as  a 
function  of  frequency  by  a  cubic  spline.  LF  takes  on  the  integer  values  1 
through  4. 

SUBROUTINE  FUNCVF (MD ,  LF) 

This  subroutine  called  by  FUNSPL  selects  for  LF  =  1(2),  the  real  (imaginary) 
part  of  the  excitation  factor  for  fitting  to  a  cubic  spline.  If  LF  =  3(4), 
the  real  (imaginary)  part  of  the  eigenangle  is  selected  for  fitting  to  a  cubic 
spline. 
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SUBROUTINE  SPLINE (X.  Y,  B,  C,  D,  N) 


This  subroutine  called  from  FUNSPL  determines  the  coefficients  B,  C,  D  of  a 
cubic  spline  interpolating  the  given  curve  [X(I),  Y(I),  I  =  1,  2,  .  .  .  ,  N] . 
If  X(I) .LE.XX.LE.X(I+1)  and  H  =  XX  -  X(I),  then  the  interpolated  value  at  XX 
is  F(XX)  =  Y(I)  +  B(I)  *  H  +  C(I)  *  H  **  2  +  D(I)  *  H  **  3.  The  interpolated 
value  is  evaluated  by  using  the  function  SPEVAL  discussed  subsequently. 

SUBROUTINE  TPLOT(FREQ.  FL,  FO,  FC,  DELTAF ,  PULSEW,  PULSED,  NUMPLS, 

IFLGTR ,  NRPT1 ,  NF) 


TPLOT,  called  from  MAIN,  sets  up  the  arrays  for  plotting  the  transmitter 
spectrum.  The  arguments  of  TPLOT  are  as  follows: 

FREQ  -  input  frequency  in  Hz. 

FL  -  lowest  frequency  in  Fourier  integral  evaluation. 

FO  -  carrier  frequency  in  Hz. 

FC  -  chip  frequency  in  Hz. 

DELTAF  -  frequency  interval  in  Hz,  at  which  transmitter  spectrum  is 
calculated. 

PULSEW  -  pulse  width  in  |js  when  used  with  IFLGTR  =  1;  and  1/e  pulse  half 
width  in  ps  when  used  with  IFLGTR  =  2. 

PULSED  -  pulse  delay  in  ps  (used  with  IFLGTR  =  1  or  2). 

NUMPLS  -  number  of  pulses  (used  with  IFLGTR  =  1  or  2). 

IFLGTR  -  transmitter  flag.  IFLGTR  =  1  corresponds  to  sequence  of  square 

wave  pulses.  IFLGTR  =  2  corresponds  to  sequence  of  Gaussian 
pulses.  IFLGTR  =  3  corresponds  to  a  differential  phase-encoded 


signal  waveform  composed  by  frequency  shift  keying  (FSK)  a 
carrier  frequency  with  a  pseudonoise  sequence  of  pulses  or 
chips.  The  FSK  index  is  0.5.  IFLGTR  =  4  corresponds  to 
Williams'  source  for  slow  wave  tail  calculation. 

NRPT1  -  number  of  frequency  points  between  FU  and  FL  used  in  Fourier 
evaluation  (FU  is  the  highest  frequency  in  Fourier  integral 
evaluations) . 

NF  -  the  number  of  frequencies  read  in. 

SUBROUTINE  TKXMTR(K,  F,  F0,  FC,  PULSEW,  PULSED,  NUMPLS,  IFLGTR,  LABELT, 

XMTR) 

TRXMTR ,  called  from  TPLOT,  evaluates  the  transmitter  spectrum.  The  arguments 
of  TRXMTR  are: 

K  -  integer  index  of  frequencies  for  which  transmitter  spectrum  is 

evaluated. 

F  -  frequency  in  Hz. 

LABELT  -  transmitter  label. 

XMTR  -  transmitter  evaluation. 

The  remaining  arguments  are  defined  in  TPLOT. 

SUBROUTINE  RPL0T(FREQ.  FL,  FI,  F2,  F3,  DELTAF ,  P,  Q,  NRPT1 ,  NF) 

RPLOT,  called  from  MAIN,  sets  up  the  arrays  for  plotting  the  receiver 
spectrum.  The  arguments  of  RPLOT  are: 

FI,  F2,  F3  -  frequencies  f^,  f^  and  f^,  respectively,  in  equation  6. 

P,  Q  -  exponents  P  and  Q  in  equation  6. 
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All  other  arguments  are  the  same  as  TPLOT.  RPLOT  calls  RECVR  and,  like  TPLOT, 
calls  the  controlling  plot  routine  PLSPEC,  which,  in  turn,  calls  BORDER, 
SYMBOL,  CURVE,  and  PLOT. 

SUBROUTINE  RECVR(K,  F,  FO,  FI,  F2 ,  F3,  LABELR,  P,  Q,  RCVR) 

Called  from  RPLOT,  this  subroutine  calculates  the  receiver  spectrum  according 
to  equation  6. 

K  -  integer  index  of  frequency  for  which  receiver  spectrum  is 

evaluated . 

F  -  frequency  in  Hz. 

FO  -  carrier  frequency  in  Hz. 

FI,  F2 , 

F3  -  frequencies  f ^ ,  f 2  and  f 3  of  equation  6. 

LABELR  -  receiver  label. 

P,  Q  -  exponents  of  equation  6. 

RCVR  -  receiver  evaluation. 

SUBROUTINE  CPPLOT(FREQ,  FL,  FO,  FC ,  DELTAF ,  NRPT1 ,  NF,  BANDW,  RHO) 

CPPLOT,  called  from  MAIN,  sets  up  the  arrays  for  plotting  the  channel  spectra 
and  the  product  spectra  consisting  of  the  transmitter,  receiver,  and  channel. 
All  arguments,  except  RHO,  have  been  previously  defined  in  TPLOT  and  RPLOT. 
The  argument  RHO  is: 

RHO  -  range  (km). 

CPPLOT  calls  CHANEL  and,  like  TPLOT,  calls  the  controlling  plot  routine 
PLSPEC,  which,  in  turn,  calls  BORDER,  SYMBOL,  CURVE,  and  PLOT  as  well  as  a 
utility  labeling  subroutine  called  PLLABL. 
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SUBROUTINE  CHANEL (F,  RHO,  CHNL) 

CHANEL,  called  from  CPPLOT,  evaluates  the  channel  spectra  for  frequency  F  (Hz) 
and  range  RHO  (km)  according  to  equation  46.  The  channel  evaluation  is  placed 
in  CHNL. 

FUNCTION  SPEVAL (XV AL .  X,  Y,  B,  C,  D,  N,  INIT) 

SPEVAL,  called  from  CHANEL,  evaluates  the  interpolating  cubic  spline  for  the 
data  [X(I),  Y(I) ,  1=1,  .  .  .,  N  at  XVAL] .  INIT  is  an  estimate  of  the 
interval  where  XVAL  lies,  X(INIT)  .LE.  XVAL  .LE.  X(INIT  +  1),  but  need  not  be 
used.  Set  INIT  =  0  if  there  is  no  estimate.  On  return,  INIT  will  contain  the 
interval  number. 


SUBROUTINE  NLOGN(N.  X,  Y,  SIGNT,  A,  B) 


Called  from  MAIN  when  INTFLG  . EQ.  0,  NLOGN  calculates  integrals  of  the  form  (s 
=  SIGNT). 


f  D 

exp[-i27CsAxl  1  |x(f )+iy(f ))  exp(i2nsf x)df 

J  A 

/•B-A 

=  I  (x(f  +  A)+iy(f  +  A))  exp (i2nsfx)df 
J  o 


(48) 


by  the  fast  Fourier  transform  technique  of  Cooley  and  Tukey  (ref  10).  This 
makes  use  of  digital  evaluations  at  the  frequencies 


f(L)  =  -  A)  ;  L  =  l,  2, 

<>n 


(49  > 


and  the  method  yields  evaluations  for  the  times 
=  B-A  ’  K  =  lf  2*  ••••»  2  * 


(50) 


Real  and  imaginary  parts  of  the  integral  are  then  stored  in  X(K)  and  Y(K), 
respectively.  The  weight  factors  and  endpoint  corrections  supplied  in  the 
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earlier  work  (ref  1)  were  in  error  and  have  been  abandoned  in  the  present 
work.  Though  the  consequent  error  in  the  output  parameter  range  of  interest 
(TAUO  to  TAUMAX)  in  reference  1  is  quite  small,  it  is  strongly  recommended 
that  the  present  program  be  used  instead  of  reference  1  even  though  the 
interest  may  be  solely  in  vertical  E  at  the  ground  produced  by  a  ground-based 
transmitter. 

The  quantity  s  =  SIGNT  takes  on  the  values  ±1  and  simply  allows  for  plus 
or  minus  transforms  js  desired.  It  should  be  observed  that  although  the 
region  of  significance  of  the  integrand  of  equation  48  may  be  quite  small,  the 
integration  limits  A,  B  may  of  necessity  be  quite  large  in  order  to  achieve  a 
desired  time  resolution  (see  eq  SO).  N  must  be  chosen  to  give  small  enough 
step  sizes  in  the  region  where  the  integrand  is  significant.  Specifically, 
step  sizes  must  be  small  compared  with  distances  (in  frequency  units)  over 
which  the  integrand  changes  appreciably. 

SUBROUTINE  FILON(N.  X,  Y,  TAU,  FU,  FL,  SUM,  SUMP) 

Called  from  MAIN  when  INTFLG  .NE.  0,  FILON  calculates  integrals  of  the  form 

(FU 

(x(f)  +  iy(f ))  exp(2nitf )df  (51) 

FL 

N 

by  Filon's  method  (ref  4),  using  2+1  point  evaluations  in  the  range 

FU  -  FL.  In  particular,  SUM  represents  the  N  point  evaluation  and  SUMP 
N*  1 

represents  the  2  +1  point  evaluation.  The  relative  error  |SUM-SUMP|/ |SUM| 

is  calculated  and  printed  in  MAIN.  As  discussed  in  the  introduction,  the 
evaluation  is  considerably  slower  than  can  be  achieved  by  using  NLOGN,  but  the 
method  is  used  to  give  a  measure  of  the  error  in  the  evaluation  as  noted 
above.  Note  that  FL  and  FU  are  chosen  on  the  basis  that  the  integrand  be 
sufficiently  small  at  the  end  points,  unlike  the  situation  with  NLOGN,  where 
the  limits  are  chosen  on  the  basis  of  the  time  resolution  required. 
Accordingly,  the  number  of  point  evaluations  (determined  by  N)  using  the  Filon 
method  can  be  appreciably  smaller  than  the  point  evaluations  using  NLOGN. 
Despite  this  tradeoff,  the  NLOGN  routine  appears  to  be  appreciably  faster. 
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SUBROUTINE  INITWF(IFLGTR,  PULSEW,  PULSED,  NUMPLS,  TAUO,  TAUMAX,  PLOTX3 
PLOTY3) 

Called  from  MAIN,  this  subroutine  calculates  the  input  waveforms  for  IFLGTR 
(eq  1  or  2).  PULSEW,  PULSED,  NUMPLS  and  IFLGTR  have  been  previously  defined 
(see,  e.g.,  TPLOT) .  TAUO  and  TAUMAX  are  the  minima  and  maxima  abscissa 
times  in  ms  on  the  output  waveform  plots.  PL0TX3  and  PL0TY3  contain  the 
output  plot  array  data. 


A  chart  showing  the  essential  flow  of  the  pulse  program  appears  below. 


DRIVER  PROGRAM 
CONTROLS  INPUT  OUTPUT 
PROCESSES  MODE  DATA. 
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V.  DESCRIPTION  OF  OUTPUT 


The  sample  output  shown  below  begins  with  the  namelist  output  followed  by 
the  run  identification.  The  mode  data  come  next.  For  each  frequency  (given 
in  increasing  order),  the  number  of  modes,  the  real  and  imaginary  parts  of  the 
eigenangle,  and  their  modal  equivalents  of  attenuation  rate  and  phase  velocity 
normalized  to  free-space  velocity  are  listed. 

The  calculated  output  of  the  pulse  shape  program  is  presented  following 

the  sample  output.  The  transmitter,  receiver,  channel,  and  product  (XMTR  * 

RCVR  *  CHNL)  spectra  are  given  as  a  function  of  frequency.  Not  all  2049 
NFFT 

(i.e.,  2  +1  with  NFFT  =  11)  lines  are  listed.  The  printout  is  controlled 

by  the  namelist  variables  NPRNT  and  INTPRT. 

Following  the  spectra  output,  output  is  presented  pertaining  to  the 
envelope  of  the  time  signature  of  the  output  waveform,  G(x,  z  ,  t;  z  ),  given 
by  equation  1.  Specifically,  the  time  TAU  in  seconds  is  given,  along  with  the 
envelope  in  dB/pV/m/kW  in  the  spirit  of  the  normalization  discussed  in  section 
II. 


The  signal  envelope  in  absolute  units  for  times  less  than  -0.3  ms  should 
be  identically  zero  (or  -°°  on  a  dB  scale) .  The  fact  that  this  is  not  so  is 
believed  to  be  associated  with  frequency  truncation  effects  and/or 
discontinuities  in  the  third  and  higher  derivatives  of  the  interpolated  mode 
data.  The  principal  outputs  of  the  program  are  the  plots.  In  the  present 
instance,  figures  1  through  6  are  generated.  These  are 

1.  Transmitter  spectrum  versus  frequency. 

2.  Receiver  spectrum  versus  frequency. 

3.  Channel  spectrum  versus  frequency. 

4.  Product  (transmitter  *  receiver  *  channel)  spectrum  versus 
frequency. 

5.  Input  and  output  envelopes  normalized  to  unity. 

6.  Waveform  output,  including  carrier  frequency  normalized  to 
unity. 
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Shown  in  figure  5  is  the  maximum  of  the  envelope  (SIGNAL  MAX)  expressed  in 
units  of  dB/|JV/m/kW.  The  right-hand  scale  on  the  plot  gives  the  envelope  in 
dB  relative  to  SIGNAL  MAX. 

Figures  1  through  6  just  discussed  apply  to  IFLGTR  =  1  (i.e.,  square-wave 
envelopes).  Figure  7  is  the  envelope  output  for  IFLGTR  =  2  (i.e.,  Gaussian 
envelope  output).  Figures  8  and  9  show  output  for  IFLGTR  =  3  (i.e.,  the 
correlator  output  for  an  MSK  format).  Figure  8  shows  the  amplxtude  of  the 
correlator  output  and  figure  9  the  phase  of  the  correlation  vector  given  by 
equation  23.  It  can  be  observed  from  the  labels  that  figures  8  and  9  apply  to 
a  different  ionospheric  profile  (i.e.,  channel)  than  that  used  for  figures  1 
through  7 . 

Corresponding  to  IFLGTR  =  4,  figure  10  shows  the  transmitter  spectrum  for 
Williams'  source,  figure  11  shows  a  representative  receiver  spectrum  used 
for  slow  tail  measurements,  figure  12  shows  the  channel  spectrum,  and  figure 
13  gives  the  product  spectrum.  Lastly,  the  slow  wave  tail  output  for  the 
Williams'  source  is  shown  in  figure  14. 
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NUM8ER  OF  PULSES  =  2 

PULSE  WIDTH  =  600-0  MICRO-SEC 

PULSE  OELAY  =  600-0  MICRO-SEC 


Figure  1 .  Transmitter  spectrum  versus  frequency. 
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Figure  4.  Product  spectrum  versus  frequency. 
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Figure  5.  Input  and  output  waveforms  normalized  to  unity. 
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Figure  6.  Waveform  output,  including  carrier  frequency. 
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Figure  7.  Envelope  output  for  IFLGTR=2. 
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Figure  8.  Envelope  output  for  IFLGTR=3. 
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Figure  9.  Correlator  phase  output  for  IFLGTR=3. 
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Figure  11.  Receiver  spectrum  used  for  slow-tail  measurements. 
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Figure  12.  Channel  spectrum. 
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Figure  13.  Product  spectrum. 
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NAVY  DEPARTMENT 
ARLINGTON.  VA  22217 
CODE  402 
CODE  420 
CODE  421 
CODE  461 
CODE  464 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CENTER 
4301  SUITLAND  RD  BLDG  5 
WASHINGTON,  DC  20390 


OFFICER-IN-CHARGE 
WHITE  OAK  LABORATORY 
NAVAL  SURFACE  WEAPONS  CENTER 
SILVER  SPRING,  MD  20910 

CODE  WA501  NAVY  NUC  PRGMS  OFF 
CODE  WX21  TECH  LIBRARY 

COMMANDER 

NAVAL  TELECOMMUNICATIONS  COMMAND 
NAVTELCOM  HEADQUARTERS 
4401  MASSACHUSETTS  AVE,  NW 
WASHINGTON.  DC  20390 
CODE  24C 

COMMANDING  OFFICER 
NAVY  UNDERWATER  SOUND  LABORATORY 
FORT  TRUMBULL 
NEW  LONDON.  CT  06320 
PETER  BANNISTER 
DA  MILLER 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
NAVY  DEPARTMENT 
WASHINGTON.  DC  20376 
NSP-2141 


COMMANDER 

ADC/DC 

ENTAFB.CO  80912 
DC  (MR  LONG) 

COMMANDER 
ADCOM/XPD 
ENTAFB.CO  80912 
XPQDQ 
XP 

AF  GEOPHYSICS  LABORATORY,  AFSC 
HASCOMAFB,  MA  01731 
CRU  (S  HOROWITZ) 

AF  WEAPONS  LABORATORY,  AFSC 
KIRTLAND  AFB,  NM  87117 
SUL  (2) 

DYC 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13440 
EMTLD  DOC  LIBRARY 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 
EEP  JOHN  RASMUSSEN 


COMMANDER  IN  CHIEF 
STRATEGIC  AIR  COMMAND 
OFFUTT  AFB.NB  68113 
NRT 
XPFS 

LAWRENCE  LIVERMORE  NATIONAL 
LABORATORY 
PO  BOX  808 

LIVERMORE.  CA  94550 
TECH  INFO  DEPT  L-3 

LOS  ALAMOS  NATIONAL  SCIENTIFIC 
LABORATORY 
PO  BOX  1663 
LOS  ALAMOS.  NM  87545 

DOC  CON  FOR  T  F  TASCHEK 
DOC  CON  FOR  D  R  WESTERVELT 
DOC  CON  FOR  P  W  KEATON 

SANDIA  LABORATORY 
LIVERMORE  NATIONAL  LABORATORY 
PO  BOX  969 

LIVERMORE.  CA  94550 

DOC  CON  FOR  B  E  MURPHEY 
DOC  CON  FOR  T  B  COOK  ORG  8000 

SANDIA  NATIONAL  LABORATORY 
PO  BOX  5800 

ALBUQUERQUE.  NM  87115 

DOC  CON  FOR  SPACE  PROJ  DIV 
DOC  CON  FOR  A  D  THORNBROUGH 
DOC  CON  FOR  3141  SANDIA  RPT  COLL 


DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON.  DC  20234 
RAYMOND  T  MOORE 

DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TELECOMMUNICATIONS 
INSTITUTE  FOR  TELECOM  SCIENCE 
80ULDER,  CO  80302 
WILLIAM  F  UTLAUT 
L  A  BERRY 
A  GLENN  JEAN 


AEROSPACE  CORPORATION 

PO  BOX  92957 

LOS  ANGELES,  CA  90009 


IRVING  MGARFUNKEL 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  RD 
BURLINGTON,  MA  01803 
RADIO  SCIENCES 


THE  BOEING  COMPANY 
PO  BOX  3707 
SEATTLE,  WA  98124 
GLENN  A  HALL 
J  F  KENNEY 

ESL,  INC 
495  JAVA  DRIVE 
SUNNYVALE.  CA  94086 
JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
PO  BOX  8555 

PHILADELPHIA,  PA  19101 

SPACE  SCIENCE  LAB  (MH  BORTNER) 

KAMAN  TEMPO 
816  STATE  STREET 
PO  DRAWER  QQ 
SANTA  BARBARA.  CA  93102 
B  GAMBILL 
DASIAC 

WARREN  S  KNAPP 

GTE  SYLVANIA,  INC 
ELECTRONICS  SYSTEMS  GRP 
EASTERN  DIVISION 
77  A  STREET 
NEEDHAM,  MA  02194 
MARSHAL  CROSS 

ITT  RESEARCH  INSTITUTE 
10  WEST  35TH  STREET 
CHICAGO,  IL  60616 

TECHNICAL  LIBRARY 

UNIVERSITY  OF  ILLINOIS 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

URBANA,  IL  61803 

AERONOMY  LABORATORY 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
J  NEWLAND 
PT  KOMISKE 

LOCKHEED  MISSILE  &  SPACE  CO,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 
E  E  GAINES 
WLIMHOF  D/52-12 
JB  REAGAN  D652-12 
R  G  JOHNSON  D/52-12 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 
PO  BOX  73 

LEXINGTON,  MA  02173 
0  M  TOWLE 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
R  HENDRICK 
F  FAJEN 

MITRE  CORPORATION 
PO  BOX  209 
BEDFORD,  MA  01730 
G  HARDING 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD 
SANTA  MONICA,  CA  90404 
EC  FIELD,  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERIC  RESEARCH  LABORATORY 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
IONOSPHERIC  RSCH  LAB 

RAD  ASSOCIATES 
PO  BOX  9695 

MARINA  DEL  REY.CA  90291 
FORREST  GILMORE 
WILLIAM  JKARZAS 
PHYLLIS  GREIFINGER 
CARL  GREIFINGER 
H  AORY 

BRYAN  GABBARD 
R  P  TURCO 
SAUL  ALTSCHULER 


RAND  CORPORATION 
1700  MAIN  STREET 
SANTA  MONICA.  CA  90406 
TECHNICAL  LIBRARY 
CULLEN  CRAIN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
DONALD  NEILSON 
GEORGE  CARPENTER 
W  G  CHETNUT 
J  R  PETERSON 
GARY  PRICE 

STANFORD  UNIVERSITY 
RADIO  SCIENCE  LABORATORY 
STANFORD.  CA  94305 
R  A  HELLIWELL 
FRASER  SMITH 
J  KATSUFRAKIS 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 
JET  PROPULSION  LABORATORY 
4800  OAK  GROVE  DRIVE 
PASADENA.  CA  91103 
ERNEST  K  SMITH 
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